%mu=1 lax-wendroff
clear
diatah=17;
a=(15-diatah):0.05:(25+diatah);
u=f0(a);
n=(10+2*diatah)/0.05+1;
for i=1:20*17
    n=n-2;
    u_temp=zeros(1,n);
    for j=1:n
       u_temp(j)= (-0.5+0.5)*u(j+2)+(0.5+0.5)*u(j)+(1-1)*u(j+1);
    end
    u=u_temp;
end
u_temp=u;
figure(1);

set(gca,'XTick',[15:1:25]);   
set(gca,'yTick',[-0.4:0.2:1.2]); 
axis([15 25 -0.4 1.2]);
hold on;
k=15:0.05:(25);
z=f17(k);
fp1=plot(k,u_temp);
fp1.LineStyle='-';
fp1.Color='r';
fp1.Marker='.';
grid on;
hold on;
fp0=plot(k,z);
fp0.Color = 'b';
string='LaxWendroff-mu=1.0.png';
title(string);
saveas(1,string);
clf(1);


